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Abstract 

In this paper, we demonstrate an elementary method for constructing new solu¬ 
tions to Bochner’s problem for matrix differential operators from known solutions. 
We then describe a large family of solutions to Bochner’s problem, obtained from 
classical solutions, which include several examples known from the literature. By 
virtue of the method of construction, we show how one may explicitly identify a 
generating function for the associated sequence of monic ic-orthogonal matrix poly¬ 
nomials {p{x,n)}, as well as the associated algebra D{w) of all matrix differential 
operators for which the {p{x,n)} are eigenfunctions. We also include some general 
results on the structure of the algebra D{w). 


1 Introduction 

A weight function is a nonzero measureable function re : M —>■ [0, oo) satisfying the 
condition that the moments J w{x)x"'dx are all finite. A weight function gives rise to an 
inner product on the space of polynomials in x, defined by 

{fix),g{x))w ■■= j f{x)w{x)g{x)*dx. 

By Gram-Schmidt, we may obtain a sequence of pairwise orthogonal polynomials {p{x, n)} 
such that p{x, n) has degree n for each integer n > 0. Such a sequence is called a sequence 
of orthogonal polynomials for w, or if w is implied, simply a sequence of orthogonal poly¬ 
nomials. 

Bochner’s problem, introduced in the paper [^, is to determine for which weights w 
the associated polynomials are a family of eigenfunctions of some second order differential 
operator. Lucky for us, Bochner provides a solution to his problem in the very same 
paper. Up to affine changes of coordinates, the only orthogonal polynomials satisfying 
this property are the classical families of Hermite, Laguerre, and Jacobi. 

The question Bochner addresses generalizes naturally to operators of higher order. 
Consider a weight function w and a sequence of orthogonal polynomials {p(x,n)} . The 
set D{w) of all differential operators for which the sequence is a family of eigenfuntions 
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forms an algebra. Moreover, this algebra is independent of the choice of sequence of 
orthogonal polynomials for w, since after all the value of p{x, n) is unique up to a scalar 
multiple for each n. In terms of D{w), Bochner’s result tells us for which w the algebra 
D[w) contains an operator of order two. However, the classification of weight matrices 
w for which the algebra D{w) contains an operator of higher order seems to be very 
hard. The classification for operators of order 4 was done by Krall 19|, but for higher 
order is still open. 

Bochner’s problem and the notion of orthogonal polynomials also extends naturally 
to matrix-valued polynomials. A weight matrix is a matrix-valued function w{x) : M ^ 
M]\f(C) which is sufficiently nice (explained below) so as to induce a nondegenerate 
matrix-valued inner product on the space of matrix-valued polynomials defined by 


{f{x),g{x))w ■■= j f{x)w{x)g{x)*dx. 


Here g{x)* refers to the Hermitian-conjugate of g{x). Again by a process identical 
to Gram-Schmidt, but with matrix values, we may obtain a sequence of polynomials 
{p{x, n)} such that p{x, n) is degree n with nonsingular leading coefficient for all integers 
n > 0 and such that {p{x,n),p{x,m))w = 0/ for m n. A sequence of matrix-valued 
polynomials satisfying these properties is called a sequence of orthogonal polynomials 
for w. 

To extend Bochner’s problem to matrix orthogonal polynomials, one should consider 
the algebra of matrix-valued differential operators, eg. operators of the form 

5 = oo(x) -|- dai{x) -!-••• + d^a£{x), 


where the ai(x) € M 7 v(C[[x]]) and 6 acts on MAr(C[[®]]) by 

f(x} ■ 5 := f{x)ao{x) + f{x)ai{x) H-h f'^^\x)ai{x). 


In particular, our differential operators act on the right. This is required because of the 
noncommutativity of the coefficients: the differential operators must act on the right in 
order to be compatible with the matrix-valued inner product j^. A function /(x) is called 
an eigenfunction of 6 if there exists a matrix A G MAr(C) such that /(x) • 6 = A/(x). 
With this notion of matrix differential operators, we may state the matrix version of 
Bochner’s problem. 

Problem 1.0.1 (Bochner’s Problem for Matrix Differential Operators). Let w{x) be 
a weight matrix, and let {p{x,n)} be a sequence of orthogonal matrix polynomials for 
w{x). When does there exists a matrix differential operator of order two for which p(x, n) 
is an eigenfunction for every integer n > 0? 

Bochner’s problem for matrix differential operators is considered in numerous papers, 
including many of the papers mentioned in the references below. For a helpful survey, see 
0- Unlike the scalar case, general classification results for Bochner’s problem remain 
elusive, even for 2x2 matrices. Many papers therefore have focused instead on providing 
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new examples of Bochner pairs {w, 5), eg. a weight matrix w and a second order operator 
6 G D{w). 

More recent papers have explored the structure of the algebra D{w) of all matrix 
differential operators for which the p{x, n) are eigenfunctions. In particular, 0 0 0 
provide examples of generators and relations for D{w) for various values of w. Even 
more examples are 0]0]j^[0j^. These examples demonstrate that the structure 
of D{w) can be nuanced and interesting, unlike in the scalar case. However, despite an 
ever increasing list of examples, more general results regarding the structure of D{w) 
remain a mystery. In particular, current methods of finding generators and relations for 
D{w) are often ad hoc, or based on computational evidence, and can involve extended 
calculation. 

The purpose of the current paper is to demonstrate how under sufficiently nice con¬ 
ditions one may use Darboux transformations to create a new Bochner pair (w, 6) from 
a known Bochner pair {w, 5). Moreover, we show that when (w, (5) arises from (w, 6) by 
a Darboux transformation then the algebras D{w) and D{w) are closely related, so that 
knowledge of the structure of D{w) leads to knowledge of the structure of D{w). As a 
result, this paper leads both to new families of examples of Bochner pairs not currently 
in the literature, as well as very efficient derivations of the structure of the algebra D[w) 
of matrix differential operators for several examples of weights w already featured in the 
literature. Lastly, we prove a couple general results regarding the structure of the alge¬ 
bra D{w), not currently featured in the literature. In particular, we prove that D{w) 
must necessarily be finitely generated over its center, and that its center must be an 
affine curve. We also prove that for N = 2, the center must be rational. 

A Darboux transformation of a differential operator <5 is a new differential operator 
5 = [iv obtained by means of a factorization of the original operator 5 = vp. For 
a short survey of Darboux transformations, see 2^, and for a great read about the 
role of Darboux transformations in the context of orthogonal polynomials satisfying 
differential equations, see 0. Further background on the Darboux transformations 


relevant to this paper can be found in 1^ 1^. Finally, for a very recent article exploring 
Darboux transformations in a noncommutative context, applying in particular to matrix 
differential operators, see 1^. This notion works equally as well for matrix differential 
operators as it does for ordinary differential operators. However, given an arbitrary 
factorization of <5, there’s no reason to expect that the transformed operator 5 will be 
in D(w) for some new weight matrix w. In order to guarantee this, we must be more 
methodical about our choice of the factorization of 5. 

The method presented in this paper relies on the application of two different adjoints 
* and I on the algebra of differential operators. The first adjoint * is the standard notion 
of the adjoint, which extends the Hermitian-conjugate on matrix-valued functions and 
satisfies 

^gmy ^ 

The second f is the formal rc-adjoint of 6, and is defined in terms of the first by 


:= w{x)6*w ^(x). 
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in some neighborhood of x = 0. The requisite properties of w are discussed further in the 
section on adjoints later in the paper. For clarity, we provide a formula for the formal 
tc-adjoint of a first order differential operator in the next example. 

Example 1.0.2. Consider the first-order matrix differential operator 6 = ao{x)+dai{x). 
Then the standard adjoint of 5 is given by 

5* = ao(x)* — a^i(x)* — dai{x)*, 

and the formal rc-adjoint of 6 is given by 

(5^ = rt;(x)[ao(x)* — a[{x)* — rc(x)“^t(;'(x)ai(x)*]tc(x)“^ — dw{x)ai{x)*w{x)~^. 

Using the two adjoint definitions we have introduced, we are now able to state 
our Main Theorem, which provides a way of obtaining a Darboux transformation of 
a Bochner pair {w,5) to a new Bochner pair {w,5). 

Theorem 1.0.3 (Main Theorem). Let {w,5) be a Boehner pair supported on (xo,xi) 
with + l)tc(x) = 0 for i = 0,1 and n > 0 and 

S = ao(x) + dai(x) + d^a 2 (x)I, 

with ao(x),ai(x) matrix-valued funetions satisfying ai(x)^ = ai(x), with a 2 (x) a sealar- 
valued funetion which is nonzero on (xo,xi), and with 5 degree-preserving. Suppose that 
there exist matrix polynomials xi(x),xo(x) satisfying the following properties: 

• deg(vi(x)) = i, det(vi(x)) ^ 0 on (xo,xi) for i = 0,1 and = vovf^ 

• tr((i;i(x)a2(x)xi(x)^)“^) G Lp‘{tx{ w{x))dx) for i = 0,1 

• (uo(x)i;i(x)"^)^a 2 (x) + (■i;o(x)i;i(x)“^)'a 2 (x) -h (xo(x)i;i(x)“^)ai(x) + ao(x) = 0 
Then the following is true 

(a) there exists a smooth, matrix-valued function f on (xo,xi) satisfying vif{vif)^ = 02 

(b) w = fwf* is a weight matrix 

(c) V ■.= dvi{x) — xo(x) and p, := —/(x)/'^(x)i^^ are degree-preserving and 6 = vp. 

(d) p{x,n) := p{x,n) ■ v is a sequenee of orthogonal matrix polynomials for w 

(e) {w,5) is a Bochner pair for 5 := piv 

(f) 

D{w) D {u~^D{w)u) n M 7 v(fl) = {u~^r]u : p G D{w), ker(z^) • rj C ker(z/)}. 
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The Main Theorem is not as general as possible. The idea is to decompose 5 as 
6 = such that v is degree-preserving and rc-adjointable. The Main Theorem 

simply provides a specific context where this holds. Even more generally, we could find 
V such that 5u = u5 for some differential operator <5, eg. a Darboux conjugation rather 
than a Darboux transformation. It is also important to note that the function / in 
the Main Theorem is not unique - and in fact there are many such functions up to a 
choice of smooth, unitary matrix-valued function. To find such an /, we split w{x) as 
u{x)u{x)*, which we may do since w{x) is Hermitian. Then u~^v~^a 2 {v'^)~^u is also 
Hermitian, and therefore may be split as hh* for some matrix-valued function h. Then 
taking / = uhu~^ completes the construction. 

Motivated by our Main Theorem, we make the following definition 

Definition 1.0.4. Let {w,6) be a Bochner pair, and let {p{x,n)} be a sequence of 
orthogonal matrix polynomials for w. We call {w, 5) a Darboux transformation of 
{w, 6) if there exist matrix differential operators p and v such that 5 = vp, 5 = pv, and 
p(x, n) := p(x, n) ■ v defines a sequence of orthogonal matrix polynomials for w. We also 
say that z/ is a Darboux transformation from the Bochner pair (rc, 5) to the Bochner 
pair {w,6). We define the associated subalgebra D{w,i',w) of D{w) by 

D{w, u, w) := {v~^D{w)v) fl Mn{Q,) = : 6 € D{w), ker(i/) • 6 C ker(z^)}. 

In general, D{w,i',w) will be a proper subalgebra of D{w). However, we should 
expect the subalgebra to contain most of D{w). In particular, we have the following 
proposition, the proof of which is found in § 13.31 

Proposition 1.0.5. Let v he a Darboux transformation from the Bochner pair {w,5) to 
the Bochner pair (w,6). Suppose 6 = 5^02 -|- dai -|-oo with det(a 2 ) ^ 0 and that there 
exists a positive integer I such that 

dimc(D({(;, n, w)i/D{w, n, w)i-i) = dimc{D{w, n, w)i+ 2 /D{w, n, w)i+i), Vi > i. 

where here D{w,v,w)i represents the linear subspace of D{w,n,w) of operators of order 
at most i. Then D{w) is a finitely generated algebra over D{w,iy,w), and is generated 
by elements of order less than i. 

In practice, the previous proposition will allow us to determine the structure of D{w) 
from the structure of D[w). 

1.1 A Brief Reading Guide 

Our Main Theorem simultaneously doubles as a source of new examples of weight matri¬ 
ces w with nontrivial D{w), and a means of obtaining explicitly the value of the algebra 
D{w). As a demonstration of this, we include explicit examples of calculations of D[w) 
from the literature, as well as extensions to higher-dimensional analogues of these exam¬ 
ples. This paper itself attempts to present some very analytic ideas in an algebraic way, 
with the motivation of attracting the algebraically minded to the problem of studying 
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the structure of D{w) from the perspective of noncommutative algebra and algebraic 
geometry. For this reason, the author includes a more in-depth background in sections 
§0 and §E21 The reader who is already acquainted with the basic ideas and notation 
of our topic, and who would appreciate a dessert-before-dinner approach to the results, 
is encouraged to skip forward to the examples in § 01 The reader who is less acquainted 
with the theory of classical orthogonal polynomials and the current theory of their ma¬ 
trix counterparts is encouraged to peruse the background sections § 12.11 and § 12.21 Of 
course, we encourage all readers to glance at § 11.21 whenever they are uncertain about 
terminology or notation. 

1.2 Terminology and Notation 

As a basic notation used throughout the paper: 

• I the identity matrix (with size determined by context) 

• C[[x]] the algebra of power series in x 

• n the Weyl algebra, ie. the algebra of differential operators with C[x]-coefficients 

• n the algebra of differential operators with C[[x]]-coefficients 

• Mn{R) the algebra of x A^-matrices with entries in R 

• D{w) the algebra of differential operators for weight matrix w (see 12.2.5]) 

• E{w) the algebra of eigenvalue sequences for weight matrix w (see 13.1.6]) 

• ker((5) the kernel of a matrix differential operator 5 € Mi^{VL) as a linear operator 
on the C-vector space M 7 v(C[[x]]) 

• i/~^ the inverse of € M]\f{Q) as a matrix pseudo-differential operator (if it exists) 

• ^(xQ,xi){x) the indicator function of the interval (xo,xi) 

We recall the dehnition of the Weyl algebra 11 for this paper. 

Definition 1.2.1. The Weyl algebra H is defined to be the set 

n = {ao(x) -|- dai{x) -!-••• + d^anix) : n G N, ai(x),..., an{x) G C[x]} 
where the products are dehned by means of the fundamental commutation relation 

xd — dx = 1. 

Note that this is slightly at odds with the typical definition of the Weyl algebra, 
whose fundamental commutation relation is dx — xd = 1. This is because, for reasons 
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explained later, our matrix differential operators will act on functions on the right. As 
an explicit example, consider a second-order differential operator 

5 = ao{x) + dai{x) + d^a 2 {x) G hi. 

It will act on a smooth function f{x) by 

f{x) ■ 5 = f{x)ao{x) + f'{x)ai{x) + f'{x)a 2 {x). 

This right action has the consequence of reversing the usual fundamental commutation 
relation of the Weyl algebra. 

2 Background 

2.1 Classical Orthogonal Polynomials and Bochner’s Problem 

This section is intended to provide a reader with a quick recap of the basic points in the 
theory of orthogonal polynomials. We begin with the definition of a weight function. 
For simplicity, we dodge technical measure-theoretic details by working with smooth 
weights. 

Definition 2.1.1. Let —oo < a <b < oo. A weight function supported on {xq,xi) is 
a nonnegative function r : M —)■ (0, oo) which is 0 off of (xq, xi) and positive and smooth 
on (xo,xi), satisfying the condition that the moments fj^x^r(x)dx < oo for all integers 
n > 0. The interval (xo,xi) on which r is nonzero is called the support of r. 

Given a weight function r(x), we define an inner product on C[x] by 


{p{x),q{x))r '■= / p{x)r{x)q{x)dx. 

Jr 

By the process of Gram-Schmidt elimination, we can construct a sequence of pairwise- 
orthogonal polynomials p{x, 0),p{x, l),p(x, 2),... with p{x, n) a polynomial of degree n. 
If we impose the additional constraint that each p{x, i) is monic, then this sequence of 
polynomials is unique. 

Definition 2.1.2. Let r(x) be a weight function. A sequence of polynomials {p{x, n)}))TQ 
is called a sequence of orthogonal polynomials for r if p{x, n) has degree n for all 
integers n > 0 and (p(x, n),p{x, m))r = 0 for all m, n > 0 with n ^ m. 

The sequence of monic orthogonal polynomials for a weight function r(x) can in 
practice be calculated recursively. 

Proposition 2.1.3 (Three-Term Recurrence Relation). Let r(x) be a weight function. 
Then for all integers n > 0, there exist constants Sn, tn & C such that 

xp{x, n) = p{x, n -I- 1) -I- SnPix, n) + tnPix, n - 1). (1) 
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This The values of Sn and may be defined explicitly in terms of the moments of 
r(x). 

Sequences of orthogonal polynomials that are simultaneously eigenfunctions of a 
second-order differential operator arise naturally in Sturm-Liouville theory. Three fami¬ 
lies of weights whose corresponding orthogonal polynomias satisfy a second-order differ¬ 
ential operator were known classically. These weights are listed in the table in Figure [TJ 
The sequences of orthogonal polynomials for these weights are referred to as the classical 
orthogonal polynomials. 


Family 

weight 

Hermite 

Laguerre 

Jacobi 

a:^e“"’l(0,oo), b> -1 
(1 - x)“(l -h x)'’l(_m), a, 6 > -1 


Figure 1: The Classical Orthogonal Polynomial Weights 


The importance and applicability of the classical orthogonal polynomials lead Bochner 
to ask and answer the following question: 

Problem 2.1.4 (Bochner’s Problem). For which weight functions r{x) do there exists 
second-order differential operators (5 € for which the sequence of monic orthogonal 
matrix polynomials p{x, n) are simultaneously eigenfunctions? 

With this question in mind, we make the following definition. 

Definition 2.1.5. The pair of data (r, e) with r a weight function and e a differential 
operator is called a Bochner pair if the sequence of monic orthogonal polynomials for 
r are all eigenfunctions for e. 

In other words, (r, e) is a Bochner pair if and only if there exists a sequence of 
complex numbers {A„}5)Tq such that p{x, n) ■ e = XnP{x, n) for all n > 0. The next 
theorem characterizes solutions to Bochner’s problem. See |29l |. 

Theorem 2.1.6. Suppose that (r, e) is a Bochner pair, with r supported on (xo,xi) for 
—oo < a < b < oo. Without loss of generality, we may assume that e = d‘^q{x) + d{l{x) + 
q'{x)) for some polynomials l{x),q{x) of degree 1 and 2, respectively. Then the following 
are true: 

(i) the weight function r{x) satisfies the Pearson equation 

q{x)r'{x)r~^{x) = l{x), x G (xo,xi) (2) 

(a) the polynomials q{x) and l{x) satisfy the boundary condition 

lim q{x)r{x) = lim l{x)r{x) = 0. 

x^a+ x^b— 
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Conversely, any weight funetion r{x) supported on an interval (a;o,a:i) satisfying (i) and 
(a) for some polynomials l{x) and q{x) of degree 1 and 2, respectively gives a solution 
(r, e) to Bochner’s problem. 

In the case that (r, e) is a Bochner pair, we can obtain the sequence of monic or¬ 
thogonal matrix polynomias for r in a far sleeker way than the recurrence relation of 
Equation ([T]) (j^ . 

Theorem 2.1.7. Let r{x) he a weight function satisfying the assumptions of Theorem 
\2.1.6[ and let q{x) and l{x) he as in the statement of the theorem. Then the sequence of 
monic orthogonal polynomials p{x, n) may he obtained by the Rodrigues formula 

p{x,n) = Cn(g'(x)"'r(x))^"'V“^(x) (3) 


for some constants Cn. Moreover, the generating function 

is analytic in (xo,xi) x I for some open interval I of 0 and is given by 


where A = A(x, t) satisfies 


F{x,t) = -tq'{\)) \ 

r[x) 


X — X — tq{X) = 0. 


(4) 


The answer to Bochner’s question is that the classical orthogonal polynomials are 
the only solutions to Bochner’s problem. This is the content of the next theorem. 

Theorem 2.1.8 (Bochnerj^). Up to an affine change of coordinates, the only weight 
functions r{x) which are solutions of Bochner’s problem are the classical weight func¬ 
tions. 


We can extend Bochner’s problem by asking what other differential operators 5 have 
the monic orthogonal polynomials of r as common eigenfunctions. Let {p{x,n)} be the 
sequence of monic orthogonal polynomials for some weight function r. The set 


D{r) = {5 € n : p{x, n) is an eigenfunction of 5 for all n} 

forms a subalgebra of the algebra of differential operators. For (r, e) a Bochner pair, it 
turns out that the answer is very straightforward. 

Theorem 2.1.9 iMiranian 2ll| 2(il|L If (r, e) is a Bochner pair, then 


D{r) = Cn{e) = C[e]. 

Remark 2.1.10. We point out that may properly contain (^^(e). For example, 

the operator e = d^{l — x^) — dx commutes with the first-order operator 5(1 — 
but the coefficients of the latter operator are not polynomial. 


The previous theorem shows that a second-order differential operator e forming part 
of the Bochner pair (r, e) generates the algebra D{r), and is unique up to a scalar 
multiple. Values of e for various weight functions are listed in the table in Figure [2j 
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Weight 

2nd-order operator 

Hermite 

Laguerre 

Jacobi 

— 2dx 

d‘^x + d{h + 1 — x) 

9^(1 — x^) + d{b — a — [a + b + 2)x) 


Figure 2: Classical Solutions to Bochner’s Problem 


2.2 Matrix Orthogonal Polynomials and Bochner’s Problem 

We next review the basic theory of orthogonal matrix polynomials and Bochner’s prob¬ 
lem. 

Definition 2.2.1. A weight matrix w supported on an interval (xo,iCi) is defined 
to be a function re : M —>■ Mi\f(S.) satisfying the condition that w vanishes outside of 
(xo,xi), that everywhere in (xo,xi) the matrix w is entrywise-smooth, Hermitian, and 
positive-dehnite, and hnally the condition that w has finite moments: 

The interval (xo,xi) where the matrix is nonsingular is called the support of w. 

Note that we work only with “smooth” weight matrices in order to avoid more 
delicate analytic considerations. For a survey on the analytic theory of matrix orthogonal 
polynomials, see 0]. Every weight matrix w defines a matrix-valued inner product (•, ■)w 
on the set Mjv(C[x]) of all x complex matrix polynomials by 

{p,q)w-= / p{x)w{x)q{x)*dx, Vp, g G M7v(C[x]). 

Jr 

Though it will not play an important role in our paper, the matrix-valued inner product 
above begets a traditional (scalar-valued) inner product by taking trace tic{{p,q)w)- 
A Gram-Schmidt type argument with our matrix-valued inner product shows that 
for each integer n > 0 there exists a polynomial pn of degree n with nonsingular leading 
coefficient, uniquely dehned up to its leading coefficient, such that {pn,Pm)w = 0/ for all 
m ^ n. We observe that orthogonality with respect to the matrix-valued inner prodnet 
is a stronger condition than orthogonality with respect to tr((p, q, )^). 

Definition 2.2.2. We call a sequence of matrix polynomials {p{x,n)} a sequence of 
orthogonal matrix polynomials (OMP) for w if deg(p(x,n)) = n for all integers 
n > 0, with non-singular leading coefficient, and 

{PnjPm)w — 0^) Vm n. 
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A OMP {p{x,n)} for w will be called monic if it satisfies the additional constraint 
that p{x, n) has leading coefficient equal to the identity matrix I for each fixed n. Once 
again, there exists a unique sequence of monic OMP {p{x,n)} for any fixed w. 

Generalizing Favard’s theorem, Duran and Lopez-Rodriguez proved that a sequence 
of polynomials {p(x, n)} being the sequence of monic orthogonal matrix polynomials for 
some weight matrix w is equivalent to the sequence {p{x,n)} satisfying a three-term 
recurrence relation under certain regularity conditions. 

Theorem 2.2.3 (Duran and Lopez-Rodriguez j^). Suppose that {p{x,n)} is a sequence 
of monic orthogonal matrix polynomials for a weight matrix w. Then there exist se¬ 
quences of complex-valued matrices {sn} and {tn} such that 

xp{x, n) = p{x, n -|- 1) -|- Snp{x, n) -\- tnp{x, n — 1), Vn > 1 (5) 


Conversely, given reasonably nice sequences of marices {sn}, {tn}, there exists a weight 
matrix w for which the sequence of polynomials {p{x,n)} defined by Equation ([5]) is a 
sequence of monic orthogonal matrix polynomials. 


Intuitively, the sequences {sn} and {tn} determine the moments of the weight ma¬ 
trix, and vice versa. In this way, there is a sort of dictionary between weight ma¬ 
trices and sequences of polynomials satisfying a recurrence relation of the above form 

SHBiniQ. 


Matrix polynomials p{x) have a natural action by matrix differential operators, ie. 
differential operators whose coefficients are matrix polynomials. This action is on the 
right in order to make it compatible with the matrix-valued inner product defined by w 
and the corresponding three-term recursion relation with coefficients on the left. In par¬ 
ticular, for polynomials ao{x),ai{x),a 2 {x) G M]y{C[x]) the second-order matrix differen¬ 
tial operator e = ao{x) -\-dai{x) -\-d'^a 2 {x) acts on a matrix polynomial p{x) € Mn{C[x]) 
by 

p{x) ■ € = p{x)ao{x) -|- p'{x)ai{x) + p"{x)a 2 {x). 


With this in mind, Bochner’s problem was reposed by Duran 


as the following 


Problem 2.2.4 (Bochner’s Problem for Matrix Differential Operators). Determine all 
weight matrices w{x) such that there exists a second-order matrix differential operator 
5 G M]\[{Q) for which the associated sequence of monic orthogonal matrix polynomials 
are simultaneously eigenfunctions. 


For reasons that will be discussed below, when such an 5 exists it may be taken to 
be tc-symmetric. This motivates the following definition. 

Definition 2.2.5. We define a (matrix) Bochner pair to be a pair {w,5) with w a 
weight matrix and 5 G M7v(D) a rc-symmetric second-order differential operator for 
which the sequence of monic OMP for w are simultaneously eigenfunctions. 

Thus we are interested in the question of what Bochner pairs exist. More generally, 
we are interested in the question of calculating the algebra D{w) of all matrix differential 
operators having the monic OMP of w as simultaneous eigenfunctions. 
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Definition 2.2.6. Let w{x) be a weight matrix, and let p{x,n) be the associated se¬ 
quence of monic orthogonal matrix polynomials. Then the algebra of matrix differ¬ 
ential operators associated to w, denoted D(w), is the set of all 6 € Mjv(fi) such 
that p{x, n) is an eigenfunction of 5 for all n. 

The behavior in the scalar case {N = 1) leads us to ask questions such as whether 
the solutions to Bochner’s problem for matrix differential operators must satisfy a form 
of Pearson’s equation or a Rodrigues-type recurrence relation. Duran and Griinbaum 
provide a partial result in this direction in the following theorem, restated in terms 
of formal adjoints. 

Theorem 2.2.7 (Duran-Griinbaum). Let w be a weight matrix. Then D[w) contains a 
second-order differential operator if and only if there exists 

5 = d‘^a2 -T dai -|- uq, 

with Oi G MAr(C[x]) and deg(ai) < i for all i satisfying 

(a) a 2 W = wa 2 

(b) w ■ 5* = oqw for 5* the formal adjoint of 6 

(c) a 2 W and {a 2 w)' — aiw vanish as x approaches the endpoints of the support of w 

In this case, {w, 5) is a solution to Bochner’s problem for matrix differential operators. 

Note that the first two conditions translate to the statement that 6 is formally iv- 
symmetric, while the third condition implies that the formal rc-adjoint of 6 is an adjoint 
(see Def. 12. ,‘1.31) . In other words, the content of the above theorem is that a weight 
matrix w has a second-order differential operator in D{w) if and only if there exists 
a rc-symmetric second-order differential operator. The above conditions imply that w 
satisfies the non-commutative Pearson equation 

2 {a 2 wy = aiwwaf (6) 

2.3 Adjoints of Differential Operators 

The study of adjoints of differential operators can be considered a subset of the study of 
adjoints of unbounded linear operators on certain well-chosen Hilbert spaces. For this 
reason, we briefly recall the dehnition of an unbounded linear operator and its adjoint. 
For a great reference on anoints of unbounded linear operators, as well as adjoints of 
differential operators, see [5|. 

Definition 2.3.1. Let TL he a Hilbert space. An unbounded linear operator on B 

is a linear function T : T>(T) —)■ V, where P(T) is a dense subspace of TL, called the 
domain of T. The adjoint of T is an unbounded linear operator T* defined on the set 

V{T*) := {y € Tf : X —7> {Tx,y) is continuous} 

and satisfying {Tx,y) = {x,T*y) for all x G D{T),y G D{T*). 
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Note that the fact that D(T) is dense is important for the existence of T*, which 
then exists as a consequence of the Hahn-Banach theorem. 

Differential operators make great examples of unbounded linear operators. An in¬ 
teresting question is when for given differential operator the adjoint is also a differential 
operator. 

Example 2.3.2. For example, consider the first-order differential operator 5 = dai{x) + 
ao(x) with ao(x),ai(x) G C[x]. Given a weight function r(x) supported on (xo,xi), we 
may view this as an unbounded linear operator on the inner product space L^(r(x)dx), 
induced by the action of 5 on C[x]. 

Using integration by parts, one finds that for polynomials p(x), q{x): 

{p{x) ■S,q{x))r = p{xi)ai{xi)r{xi)q{xi) - p(xo)ao(xo)r(xo)g(xo) 

+ {p{x),q{x) ■ {-ai{x)*r{x)dai{x)~^ -Fao(x)))^. 

Thus if ai(xo)r(xo) = ai(xi)r(xi) = 0, we see that the adjoint of S is the differential 
operator (51' = —ai(x)*r(x)5ai(x)“^ -|- aQ(x). In particular, if r(x) is constant on its 
support, and ai(xi) = 0 for i = 0,1, then the adjoint (5^ of <5 is <5^ = —ai(x)*9 -|- ao(®)- 

A similar observation applies to matrix orthogonal polynomials, and this motivates 
the definition of a ^-operation on the algebra of matrix differential operators. 

The algebra of matrix differential operators with coefficients in the power se¬ 

ries ring C[[x]] is equipped with a canonical ^-operation, extending the usual ^-operation 
on Mj\f{C[[x]]) defined by Hermitian conjugates and satisfying 

{xl)* = xl, {diy = -dl. 

Definition 2.3.3. Let tc be a weight matrix supported on an interval (xq, xi) containing 
0 and let 5 G A formal adjoint of 6 with respect to w, or formal tc-adjoint 

of 6 is the unique differential operator (5^ G M]\f{Q,) defined on (xo,xi) by 

(5^ := w6*w~^. 

(Note that this definition relies on the assumption that w is Hermitian.) In particular, 
the ^-operation is exactly the formal tc-adjoint for w{x) = The assumption 

that 0 is in (xq, xi) may always be achieved by means of an affine change of coordinates, 
and is necessary for w to have a local representation as a power series based at 0. 
Moreover, the assumption that w is positive-definite also implies that w~^ has a power 
series expansion at 0. Consequently, the formal tc-adjoint is indeed an element of M]\[{Q). 
As a notational point, we will use * to denote the canonical adjoint on Mn{Q) and will 
use t to denote the formal adjoint with respect to a particular weight matrix w (the 
value of w will be implied from the context). 

Definition 2.3.4. Let re be a weight matrix, and <5 G Mn{^) a differential operator. 
Then 6 is called formally tc-symmetric if (5^ = (5. 


13 



We next define the adjoint of a differential operator. To do so, we must define the 
Hilbert space upon which it acts. 

Definition 2.3.5. Let w he a weight matrix. We define the Hilbert space of w to be 

^{{w) := MN{L‘^{tr{w)dx)). 

We call a matrix differential operator 5 € amenable if p{x) ■ 6 G 1-L{w) for all 

matrix-valued polynomials p{x) G Mn{C[x]). 

Definition 2.3.6. Let whe a weight matrix and let S G Mn{^) be an amenable matrix 
differential operator. An amenable matrix differential operator 6 G MAr(H) satisfying 
the identity 

{p-6,q)w = {p,q-S)w, Vp,q G Mn{C[x]), 

is called an adjoint of 5 with respect to w, or w- adjoint of 6 . If 5 = 6 , then 5 is 
called tc-symmetric. 

Every differential operator has a formal rc-adjoint, but not necessarily a tc-adjoint. 
Put another way, even though each differential operator 5 will have an adjoint as an 
unbounded operator on T-L{'w), this linear operator adjoint need not be a differential 
operator. 

Proposition 2.3.7. Let w be a weight matrix and let 6 G M]\f(Ll) be an amenable matrix 
differential operator. Then an adjoint of 5 with respeet to w, if one exists, is equal to 
the formal w-adjoint 5^. 

Example 2.3.8. Consider the weight function r(x) = e“®l(o,oo)(a^) supported on the 
interval (0, oo). The formal r-adjoint of the differential operator d is 

dt = e-^d*e^ = e-^{-d)e^ =-8+1. 


However, 

and 


(1-5,1), = (0,1), = 0 


(l,l-5t), = (1,1), = 1. 

Therefore 5^ is not an r-adjoint for 8 and thus 8 has no r-adjoint. 

As another example, the formal adjoint of the operator 8x is 

{8x)^ = e~^{8x)*e^ = e~^{—x8)e^ = —8x + x — 1. 

Moreover, integration by parts shows us that 

{p'{x)x, q{x))r = {p{x), -qfx)x + q{x){x - 1 )), 

for all polynomials p{x),q{x) G C[x]. Therefore —8x + x — 1 is in fact an r-adjoint of 8x. 
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Expanding on these examples, we have the following lemma 

Lemma 2.3.9. Let v = dfi{x) + /o(x) he a matrix differential operator for /i,/o G 
M]\f{C[x]). If fi{xi)w{xi) = 0 for i = 0,1 then v has a w-adjoint. 

Proof. For any polynomials p{x),q{x) G M]\f{C[x]), integration by parts tells us that 

{p{x) ■ u, q{x))^ = p{x)fi{x)w{x)q*{x)\%l + {p{x), q{x) ■ 

Thus if fi{xi)w{xi) = 0 for i = 0,1 the formal zn-adjoint of v is the ra-adjoint. □ 

The application of adjoints of differential operators to our situation is provided by 
the following theorem. 

Theorem 2.3.10 (Griinbaum-Tirao [I^). Let w be a weight matrix, and let 6 G D{w). 
Then a w-adjoint of S exists and is in D{w). The operator (i i—>■ <5^ is an involution on 
D{w) giving D{w) the structure of a ^-algebra. 

As a consequence, we have the following corollary originally proved in Hi- 

Corollary 2.3.10.1. Let w be a weight matrix. Then D{w) contains a differential 
operator of order m if and only if D{w) contains a w-symmetric differential operator of 
order m. 

Proof. If D{w) contains a second-order matrix differential operator ca, then by the previ¬ 
ous theorem D{w) contains uj and i{u} — ca^). Both operators are clearly symmetric, 
and by definition of the adjoint of order at most m. Since 

1 1 

UJ = -{UJ + UJ^) - -iii{uj - w'^)) 

and oj is of order m, at least one of the two symmetric operators must be order m. □ 

3 Proof of the Main Theorem 

3.1 Degree-Preserving Differential Operators 

Definition 3.1.1. We call a matrix differential operator 6 G degree-filtration 

preserving if for all polynomials q{x) G Mi\f(C[x]), the degree of q{x) • 5 is no greater 
than the degree of q{x). We call 5 degree-preserving if the degree of q{x) ■ 6 is equal to 
the degree of q{x) for all q{x) G Mv(C[x]). In particular, degree-preserving differential 
operators necessarily act injectively on the algebra of matrix polynomials. 

We use to denote the subalgebra of Mjv(fl) of degree-filtration preserving 

matrix differential operators. Due to its distinguished role in the following, we fix the 
notation s = dx. 

Lemma 3.1.2. The subalgebra is equal to the Mn{C)- subalgebra of Mjv(D) 

generated by s and d 
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Proof. Note that both s and d are degree-filtration preserving, and therefore the M]\f{C)- 
subalgebra of Mn{Q) that they generate is contained in . Thus to prove our 

lemma, it suffices to show the opposite containment. Suppose that 6 is degree-filtration 
preserving, of order n. Since 6 G we know that 

n 

S = ^ d'^ai{x) 

i=0 


for some ai{x) G M]\f{C[x]). 

We claim that deg(aj(x)) < i for all 0 < i < n. To see this, suppose otherwise. Then 
let j be the smallest nonnegative integer satisfying degaj{x) > j. Then 


X- 


■n • (5 = 


•J -I 

^ (j — 0- 

i=o ’ 




ai{x) 


is a polynomial of degree greater than j, contradicting the assumption that 5 is degree- 
filtration preserving. This proves our claim. 

Next note that for all integers j > 1, 


d^x^ = s{s — l)(s — 2)... (s — j -|- 1). 

Therefore if a{x) G Mj^{x) is of degree < j, then d^aj{x) is in the M 7 v(C)-subalgebra of 
generated by d and s. From this it follows that 6 is in the subalgebra of 
generated by s and d. This proves our lemma. □ 

Remark 3.1.3. The expression d^x^ = s(s — l)(s — 2)... (s — j — 1) found in the 
proof above is exactly the reason why the substitution x = may be used to re¬ 
express an Euler-Cauchy equation as a linear ordinary differential equation with constant 
coefficients. 


Proposition 3.1.4. The natural map 

OO 

09”MAr(C[s]) 

n=0 

is an isomorphism of Mn{C[s])- modules. In particular, every 5 G Mn{Q)^ has a natural 
expression of the form 

n 

6 = ^d"ai(s) 
i=0 

for some integer n > 0 and some matrices ai{s) G M 7 v(C[s]). 

Proof. Note that sd^ = d^{s + j), and therefore a{s)d^ = d^a{s + j) for all a(s) G 
M 7 v(C[s]). Using this relation, the above follows. □ 


16 



Proposition 3.1.5. Let w be a weight matrix, and let p{x,n) he a sequence of manic 
orthogonal matrix polynomials for w. If 5 ^ D{w), then 5 is degree-filtration preserving. 
Consequently S = XlILo some ai{s) € Mjv(C[s]), in which case 

p{x, n) ■ 6 = ao{n)p{x, n), Vn E Z+. 

Proof. Suppose that S E D{w). Then for all integers i > 0, there exists E M^iC) 
such that p(x, i) • 5 = \ip{x,i). Furthermore, since the p(x, i) form a basis for MAr(C[x]), 
given a q{x) E Mj\f{C[x]) of degree n, we may write q{x) = ^11=0 p) some 
Co,... Cn E Then 

n 

q[x) ■ 8 = '^Ci\ip{x,i), 
i=0 

which has degree at most n. Since q{x) was arbitrary, this shows that 5 is degree- 
hltration preserving. Therefore by the previous proposition, 8 = J2i=o ^^ai(s) for some 
ai(s) E MAr(C[s]). 

We know that p{x,n) = Ix'^ + (lower degree terms), and therefore 

r 

p{x,n) ■ 8 = ^^p(x,n) • 5*Oi(s) = ao(n)x"' + (lower degree terms). 
i=0 

Therefore since p{x,n) ■ 8 = Xnp{x,n) = + (lower degree terms), we have that 

An = ao{n) for all integers n > 0. □ 

The previous proposition in particular shows that if p{x, n) is a sequence of monic 
orthogonal matrix polynomials for a weight matrix w and if E D{w), then the sequence 
{An}^o — ^n{C-) satisfying Xnp{x,n) =p{x,n) ■ 8 will be a polynomial in n. Thus we 
may dehne a map A : D(w) —>■ M 7 v(C[n]) satisfying the property 

A(8)(n)p(x, n) = p{x, n) ■ 8. 

We denote the image of the map A in Mj\f{C[n]) by E{w). 

Definition 3.1.6. We call the subalgebra E{w) the algebra of eigenvalue sequences 
associated to the weight matrix w. We call the map A the eigenvalue isomor¬ 
phism. 

The eigenvalue homomorphism defines an njection of D{w) into Mv(C[n]). This was 
shown in [1^, but is reproved here in our notation. 

Proposition 3.1.7. The map A is injective, and in particular defines an isomorphism 
of D{w) onto E{w). 

Proof. Suppose that A{8) = A{8') = A(n) E M 7 v(C[n]) for some 8,8' E D{w). Then this 
in particular implies that for all n, 

p{x, n) ■ {8 — 8') = (A(n) — X{n))p{x, n) = 0. 

It follows that the kernel of 5 — (5' contains Mm{C.){pq,pi, ... } = Mj\f{C[x]). Therefore 
8 — 8' must be identically zero. □ 
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3.2 Kernels of Differential Operators 


In this section we touch on some necessary results regarding the kernel of a matrix 
differential operator of order r. By default throughout the paper, unless stated otherwise, 
the kernel of a differential operator always refers to its kernel as a linear operator on the 
ring of power series M 7 v(C[[x]]). 

Definition 3.2.1. Let fi € be a matrix differential operator. We define the 

kernel of fj, to be the subspace ker(|u) C Mjv(C[[x]]) defined by 

ker(|u) = {V’ G Mn{C[[x]]) : ip ■ fi = 0/}. 


For this reason, the elements of the kernel may not in fact belong to the Hilbert 
space determined by the weight matrix w. Furthermore, in this section the term matrix 
differential operator will always refer to an element in Mi\f{Q). 


Lemma 3.2.2. Suppose that p, is a monic matrix differential operator of order r. 
Then ker(;u) is a loeally free as a left Mp^{C)-module, generated by some ipi,... ,ipr G 
Mj\f{C[[x]]) with Ipi a unit in MAr(C[[x]]) for all i. 

Proof. There exist ao(x),... ,ar-i{x) G MAr(C[[x]]) such that p = d''I + X]i=o d’^ciiix). 
Then the matrix differential equation 


r—1 


+ '^y^"\x)ai{x) = 0 
i=0 


may be linearized to 




^ 0 

0 . 

. 0 

-ao(x) 



7 

0 . 

. 0 

-ai(x) 

Y' = YU, y= ( y y' . 

.. ), U = 

0 

7 . 

. 0 

-a2(x) 




0 . 

. 7 

-ar-i{x) 


for Y G U G C[[x]]''^''’’^): Let Uo,Ui,--- G M^jv(C) with U = Y.ZoUix\ 

Then Y = ^ solution, for all c* G satisfying the recursion relation 


1 

Cj+i = ^ CiUj-i. 

i=o 

This establishes a M 7 v(C)-linear map between the space of solutions Y to T' = YU and 
the space of constant matrices cq G , which is a free Mjv(C)-module. Under this 

isomorphism, the basis 


(7,0,0,...,0), (7,7,0,...,0), (7,0,7,...,0), (7,0,0,... , 7), 


corresponds to a MAr(C)-module basis ipo,... ,ipr of ker(^). By choice, each of the ipi 
satisfies ipi{0) = 7, and therefore ipi is a unit in M 7 v(C[[x]]) for all i. □ 
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Remark 3.2.3. A point of caution is warrented here. It is important in the above 
theorem that the coefficients of the operator reside in the power series ring C[[x]]. For 
example, the monic differential operator d — has no kernel in C[[x]], or even in the 
ring of Laurent series C((x)). 

Remark 3.2.4. The assumption that ^ is monic is also important. For example, if we 
consider the differential operator 


= 



1 

0 


5 , 


the kernel of /x is , which is not a finite-dimensional vector space. 

Lemma 3.2.5. Suppose that p, is a monic matrix differential operator of order r, and 
that ker(/x) has a Mj^{C)-module basis consisting of units Then p has a factorization of 
the form 

p = {dl - (pi){dl -(p2)-- - {dl - (fr)- 

for some (pi,... ,ipr G MAr(C[[x]]). 

Proof. We proceed by induction on r. If r = 1, then the statement is trivial. Assume 
the result for all monic operators of order r — 1, let /U be a monic operator of order r, 
and choose V' G ker(|u) such that V’ is a unit in M 7 v(C[[x]]). Writing 

r—1 

'il)p'if~^ = d''! -|- d'^bi{x) 
i=0 

for some bifx) G C[[x]], the fact that I ■ {'ijjp'ip~^) = 0 implies that bo(x) = 0. Therefore 
'ipp'ilj~^ = dv for some monic matrix differential operator v of degree r—1. Note that 
is monic of order r — 1, so by our inductive assumption, it has the desired 
factorization, and therefore 

p = vf) = {dl — 

has the desired factorization also. □ 

The previous lemma allows us to provide necessary and sufficient conditions for a 
monic matrix differential operator to divide another matrix differential operator on the 
left. 

Lemma 3.2.6. Suppose that p,S are matrix differential operators, with p monie. Then 
p~^5 is a matrix differential operator if and only ifkei{p) C ker((5). 

Proof. From the previous lemma, it suffices to consider the case when p is monic of 
order 1. In this case, ker(/x) is generated as a left iWn(C)-module by £ MAr(C[[x]]) 
with a unit in M]\f{C[x]). Note that p — {dl — is a differential operator of 
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order 0, and has kernel containing ip. Since V' is invertible in Mi\f(C[[x]]) it follows that 
= dl — 

Now since ^p G ker((5), we have that is in the kernel of . It follows 

by the same argument as in the previous lemma that = dr] for some differential 

operator rj. Therefore 

5 = 'ip~^dip'ip~^r]'tlj = {d — = ]M'ip~^ rfip. 

Therefore is a matrix differential operator. 

Conversely, suppose := is a matrix differential operator. Then for any G 
ker(z^) we have that 


-ip ■ 5 = Ip ■ {v]j) = {ip ■ u) ■ ]j, = 01 ■ ]j, = 01. 

Therefore V’ G ker((5). This completes the proof. □ 

We can loosen the restriction that ]i is monic by instead assumin gthat there exists 
ly G Mjv(fl) such that I'fi is monic. This leads to the following generalization of the 
previous lemma. 

Lemma 3.2.7. Suppose that ]i, 5 are matrix differential operators, and suppose that 
there exists a matrix differential operator v such that v]i is monic. Then ]i~^5 is a 
matrix differential operator if and only i/ker(/i) C ker((5). 

Proof. Clearly if ]j,~^5 is a matrix differential operator, then ker(;u) C ker((J). Thus it 
suffices to prove the converse. Suppose that ker(/i) C ker((5). Then ker(z/|u) C ker(z^(5), 
and therefore by the previous lemma, ]i~^5 = = ]i~^5 is a matrix differential 

operator. □ 

Instead of left-dividing, we are sometimes interested in conjugation. The next lemma 
provides necessary and sufficient conditions for the conjugate to be a differential operator. 

Lemma 3.2.8. Suppose that ]i, 5 are matrix differential operators, and suppose further¬ 
more that there exists a matrix differential operator u such that v]i is monic. Then 
is a matrix differential operator if and only f/ker(;u) • S C ker(^). 

Proof. From the previous lemma, ]i~^5]i is a differential operator if and only if ker(/r) C 
ker((5^), which in turn holds if and only if ker(^) • 6 C ker(^). □ 

Remark 3.2.9. In the Terminology and Notation section above, we talk about ]i~^, de¬ 
fined as a matrix pseudo-differential operator (if it exists). By matrix pseudo-differential 
operator, we mean a pseudo-differential operator with coefficients in Mi\f{C{{x))) for 
C((x)) the field of Laurent series in x. When does this pseudo-differential operator 
exist? The answer is that ]i~^ exists if and only if ]i is not a zero-divisor in 
In particular, this holds if the leading coefficient of ]i has a determinant which is not 
identically zero. 
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Remark 3.2.10. Suppose that we know exists and is, in fact, some differential 

operator by the above lemmas. How can we calculate its precise value? One way is to 
note that 5 := must be a differential operator of order at most (. + 2m where 

i is the order of 5 and m is the order of v. The reason for this is that the order of 
is bounded by the order of u. Thus we can write 5 = some 

unknown functions an{x) and then determine which values of the an{x) work to satisfy 
the equation Su = ud. This removes the requirement of calculating explicitly. Note 
also that if v has nonsingular leading coefficient, then the order of v~^5v is the same as 
the order of 5, and therefore in such a case we may assume an{x) = 01 for n > i. 

3.3 The Proof 

Before proving the Main Theorem, we introduce one more lemma to help with some 
moment estimates. 

Lemma 3.3.1. Suppose that a,b ^ Mjv(C) are positive-semidefinite. Then 

tr(a6) < tr(a)tr(6). 


Proof. By Cauchy-Schwartz, 

tr(a 6 ) < (tr(a2)tr(62))V2_ 

Furthermore, since a is positive-semidefinite, 

tr(a^) < tr(a)^ 

and similarly for b. From this the statement of the lemma follows immediately. □ 

We now have everything in place for the proof of the Main Theorem (Theorem 

fnm i. 

Proof of Main Theorem. 

(a) Note that since det(u(x)) € C\{0}, is also a polynomial. Moreover since 

w{x) is smooth and positive-definite on its support (xo,xi), we may factor w{x) = 
u{x)u{x)* for some smooth function u{x) of full rank on (xo,xi). Next, since 5 
is u)-symmetric, the leading coefficient 02 must be rc-symmetric also. Therefore 
is also rc-symmetric, and consequently u~^v~^a 2 {v~^)'^u is smooth and 
Hermitian on (xo,xi). Thus we may factor it as 

u~^v~^a2{v~^)^u = hh* 

for some smooth function h on (xo,xi). Taking / = uhu~^, we have that 
fp = uhu~^w{uhu~^)*w~^ = uhh*u~^ = 
and therefore / satisfies the properties stated in (a) 
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(b) From the definition of /, we have that | det(/)p = | det(/i)p = | det(?;)|“^ det(a 2 ) / 0 
on (xo,xi). Therefore f{x) has full rank for all x € (xo,a:i), and it follows that 
w{x) = f{x)w{x)f{x)* is positive-definite for all x G {xo,xi). Moreover, w{x) is 
smooth on (xo,xi), since it is the product of three smooth matrix-valued functions. 
Therefore to prove that w{x) is a weight matrix, all that is left to show is that w{x) 
has finite moments. Equivalently, we must show 



x\'^ti{'w)dx < oo. 


To show this, first note that 


tr({(;) = tr(fuu*f*) = tr{uhh*u*) = tT{hh*u*u). 

By Lemma (1,1.3. ip 

tr{hh*u*u) < tr{hh*)tr{u*u) = tr(ii“^u“^a 2 (u“^)^rt)tr(t(;) 
= tr((u“^)^a2u“^)tr(t(;) 


and therefore 

rxi 


l-Xl l-Xl 

/ \x\^tT{w)dx < / |x|”'tr(u“^(x)a 2 (x)(u“^(x))^)tr(tc(x))(ix < oo 

JXQ JXQ 


by Holder’s inequality and the fact that w{x) has finite moments. Therefore tc is a 
weight matrix. 


(c) We calculate 

= -v^^a2{d - 
= —Vi^a2{—d — w'w~^ — 

= ^02 -I- v'^^a2w'w~^ -(- uf^a2(uouf^)'*' -|- (uf^ 02 )' 


Furthermore, the noncommutative Pearson equation tells us that 

(uf^a2)' = (uf^)'a2 -I- + ®i) “ a2w'w~^. 

Substituting this in and simplifying, we find 

/i = 5 u]”^02 - Uq'^oo -h UQ'^[(uouf^a2(uouf^)^ - {vQv'{^fa2] 

Then using the fact that (uouf the last summand cancels out, leaving 

/X = dv^^a2 - Uq ^ao 
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Using this, we calculate 


ly/i = [dvi - Uo)((9Ui ^02 - Uq ^ao) 

= 5^02 + d[{vi)'v^^a2 - UQuf ^02 - uiUq ^oo] + ao 
= 9^02 + 9[-(t'ou7^)“^(uou7^)'a2 - uouf^a2 - + oq 

= 5^02 + d - (uouf^)“^[(uouf^)'a2 + (uouf^)^a2 + oq] + oq 

= 5^02 + dai + oo = (5. 

Next note that v is degree-filtration preserving. If v is not degree preserving, then 
there exists a nonzero polynomial p{x) satisfying p{x) ■ v = 0. Then 

deg(p) = deg(p • 6) = deg{p ■ u ■ = deg(O) = —oo. 

This is a contradiction, and therefore is degree-preserving. Since 6 = vp and v are 
degree-preserving, it follows that p is also degree-preserving. 

(d) The assumption that \\m.x^xi{\x\^ + l)rc(x) = 0 implies that vf is rc-adjointable by 
Lemma (I2.3.9p . Moreover, we calculate: 

(p(x,m),p(x,n ))^5 = {p{x,m) ■ vf,p{x,n) ■ vf)y, = {p{x,m),p{x,n) ■ 

= {p{x,m),p{x,n) ■ {up))w = {p{x,m),p{x,n) ■ 6)w 

Since 5 is degree-filtration preserving, if m > n then p{x, n) ■ {5 + vc) may be written 
as a MAr(C)-linear combination of p{x, 0),... p{x, n), and therefore {p{x, m),p{x, n) ■ 
{5 + i'c))w = 0. Consequently (p{x,m),p{x,n))w = 0/ for m> n. If m < n, then 

{p{x,m),p{x,n)){s = {p{x,n),p{x,m))*~ = (0/)* = 01. 

Therefore in any event 


{p{x, m),p{x, n))q, = 01, n. 

Moreover p{x, m) is a polynomial of degree m for each m, since is degree-preserving. 
It follows that p{x, m) is a sequence of orthogonal matrix polynomials for w. 

(e) We know that p{x,n) ■ 5 = \np{x,n) for some \n G Mj\f{C) for all n > 0. We 
calculate 

p{x, n) ■ 6 = p{x, n) ■ {vpv) = p{x, n) ■ {5v) = Xnp{x, n) ■ v = \np{x, n). 

Since p{x, n) is a sequence of orthogonal matrix polynomials for w, this shows that 
6 G D{w). To complete the proof, we must show that 5 is ui-symmetric. Since 
6 G D{w), we know that S is ui-adjointable, so it suffices to prove that 5 is formally 
ui-symmetric, ie. that 

5i((5)*(5i)“^ = 6. 
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We calculate 


= fwrifpiy^urirr^w-^r^ 

=/u;/V*(z.t)*(/t)*r(/*)-i^-i/-i 
= fwiil^f)\l'f))*W~^f~^ = /((i^/)^(l^/))V“^ 

= = s. 

This proves (f). 

(f) Suppose that 5' G {i'~^D{w)iy) fl MAr(O). Then there exists 5' G D{w) satisfying 
(5V = 1/5'. Moreover, there exists a sequence Aq, A'j^, • • ■ G Mn{C) satisfying p{x, n) ■ 
5' = \np{x, n) for all n > 0. Therefore we have that 

p{x, n) ■ 5' = p{x, n) ■ uS' = p{x, n) ■ {5'v) = A^p(x, n) ■ v = XnP{x, n). 

This shows that 5' G D{w), proving that {v~^D{w)i') D Mn{Q) C D{w). The fact 
that {u~^D{w)u) n Mn{Q) = {i>~^5i> : 5 G D(w), ker(i/) • 6 C ker(z^)} is just a 
restatement of Lemma (13.2.81) . 

□ 

We next prove Proposition 11.0..^ Before doing so, we establish the following lemma. 

Lemma 3.3.2. Let v be a Darboux transformation from a Bochner pair {w, 5) to a 
Bochner pair (w,5). Then D{w)5 C D{w,u,w). 

Proof. Suppose that rj G D{w). Then there exists p G M]\f{Ll) satisfying 5 = i/p and 
5 = pu and therefore 

p := i'{r]5)i/~^ = i/rjp G M]\f{Ll). 

Let {p{x,n)} be a sequence of orthogonal matrix polynomials for w. Then there exists 
a sequence Ao,Ai,--- G satisfying p{x,n) ■ 5 = Xnp{x,n) for all n > 0. Also 

p{x, n) := p{x, n) ■ v defines a sequence of orthogonal matrix polynomials for n;, and 
therefore there exists a sequence Ao,Ai,-'' G M]\f{C) satisfying Xnp{x,n) = p{x,n) ■ rj. 
We calculate 


p{x, n) ■ p = p{x, n) ■ (i/pp) = p{x, n) ■ {pp) 

= Xnp{x, n) ■ p = Xnp{x, n) ■ [vp] 

= Xnp{x,n) ■ 5 = XnXnp{x,n). 

Therefore p G D{w). Furthermore iy~^pv = p6, and it follows that p5 G D{w,i/,w). 
Since p G D{w) was arbitrary, this proves our Lemma. □ 
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Proof of Proposition \1.U. 51 Suppose that 5 and i satisfy the assumptions of the state¬ 
ment of the proposition. Then for all rj G M]\f{PL), the order of ^6 is the order of r/ plus 
two. 

By the previous lemma, D{w)5 C D{w,v,w). Therefore multiplication by 5 defines 
a C-linear monomorphism D{w)i D{w, v, rc)i+ 2 - This in turn restricts to a monomor¬ 
phism 

D{w)i/D{w)i_i D{w,v,w)i+ 2 /D{w,v,w)i+i. 

Furthermore, the inclusion D{w,i',w) C D{w) induces an injection 

D{w,i^,w)i/D{w,i^,w)i-i D{w)i/D{w)i_i 


and therefore 


dime ( 


D{w)i 


< dime 


D{w,v,w)i+2 


^D{w,i',w)i-i J \D{w)i-iJ \D{w,u,w)i+i 

Thus for i> i aW. of the above dimensions are equal, and therefore 

D{w,i',w)i/D{w,L',w)i-i A D{w)i/D{w)i_i, Vi > i. 
Hence D{w) is generated over D{w,i',w) by elements of order < i. 


□ 


4 Explicit Examples 

In this section we provide explicit examples of the Main Theorem in action. The method 
of finding examples is straightforward: we consider a specific Bochner pair (w, 5) satis¬ 
fying the assumptions of the Main Theorem, and then attempt to find a degree 1 matrix 
polynomial vi {x) and a constant matrix vq satisfying the differential equation presented 
in the Main Theorem. Once we find such a pair uo,ui(x), all that is left is to check 
that it satisfies the remaining assumptions of the Main Theorem. One helpful thing 
to point out in our search for vo,vi{x) is that if vi{x) satisfies the assumptions of the 
Main Theorem, then voVi{x)~^a 2 {x) will also be a polynomial of degree one. Therefore 
it makes sense to restrict our search to polynomials of the form 

uoui(x)"^ = {Aix -h Ao)a2ix)~^ 

such that ui(x) is also a polynomial of degree 1. Inserting this into the differential 
equation in the Main Theorem then results in a pair of algebraic relations that Aq and 
Ai must satisfy. To illustrate this point further, we provide the following two examples. 

Example 4.0.3. Consider a Bochner pair of the form (e~^^ ,6) for 6 = —d^I + d2xl + B 
for some Hermitian matrix B G Any such matrix B will still give a Bochner pair, 

so we will leave it as flexible for now. Suppose that vov{x)~^ = Aix -|- Aq satisfies the 
differential equation of the Main Theorem. Plugging this in and comparing coefficients 
of similar powers of x, we obtain the relations: 

{Ai — 2I)Ai = 0/, AjAq + AqAi — 2Aq = 01, Aq + Ai = B. 
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Choosing Aq, Ai, and B so as to satisfy these equations, we obtain a potential candidate 
for v{x). Below we explore one family of solutions of these equations. 

Example 4.0.4. Consider a bochner pair of the form (x^e“*l(o,oo))'^) ^ — —d^xl + 

d{x — (5 + 1))/ + B for some Hermitian matrix B € M]\f{C). Again, any such B 
will give a Bochner pair, so we may use B as an additional variable. Suppose that 
uou(x)“^ = {Aix + Ao)(l — x'^)~^x~^ satisfies the differential equation of the Main 
Theorem. Plugging this in and comparing coefficients of similar powers of x, we obtain 
the relations: 

Ai{Ai — I) = 01, Aq{Aq + bl) = 01, B = AiAq + AqAi + bAi — Aq. 

Choosing Aq, Ai, and B so as to satisfy these equations, we obtain a potential candidate 
for v{x). We do not explore this further in the examples. 

Example 4.0.5. Consider a Bochner pair of the form ((1 — x^Y''^,6) for 5 = —d^{l — 
+ dx{r + 2)1 + B for some Hermitian matrix B G Again, any such matrix 

B will give a Bochner pair, so we may use B as an additional variable. Suppose that 
uou(x)“^ = {Aix + Aq){1 — x^)~^ satisfies the differential equation of the Main Theorem. 
Plugging this in and comparing coefficients of similar powers of x, we obtain the relations: 

Af — r Ai + Aq = 01, AiAq + AqAj — tAq = 01, B = Ai + Aq. 

Choosing Aq, Ai, and B so as to satisfy these equations, we obtain a potential candidate 
for v{x). Below we explore one family of solutions of these equations. 

The relations in each of these examples are the central calculation which generated 
the examples provided below. Many of the results in the example section below were 
double-checked using python code with the Sympy symbolic computation library j^. 


4.1 A Family of Examples of Hermite Type 

As our first example, we consider Darboux transformations of Bochner pairs {w, 5) for 


6 = -d'^I + d2xI + B, w = e-^^I, B = 


CC* + 21 
01 



for some nonsingular, normal C G M 7 v/ 2 (C) (we need C to be normal so that 6 is 
tc-symmetric). For this to work, we must assume N is even. According to the Main 
Theorem, we can find examples by finding an invertible, degree 1 polynomial v{x) sat¬ 
isfying a certain differential equation, along with a handful of other properties. Define 
ui(x) and vq{x) by 


v{x) ^ 


21 

01 \ / 

01 

c 


x + 



01 

0/ / \ 


01 


vq{x) = I. 


Then vq{x),vi{x) satisfy the requirements stated in Theorem 1.0.3. for the Bochner pair 
{w,6). 
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Proposition 4.1.1. The functions vo{x),vi{x) as defined above satisfies the assumptions 

2 

of the Main Theorem for the Bochner pair (e“* I, (5) also defined above. 


Proof. We calculate that (ui(x) ^ — 21xvi 


CC* 01 
01 C*C 


and therefore 


i;i(x) = 


01 c*c 


(xi(x) ^—2x7) = 


0/ -2{C*C)- 


01 (c*)- 


In particular ui(x) is a polynomial of degree 1 and det(u(x)) G C\{0} (because ui(x) is 
a unit in Mn{C[x])). We also see that ui(x)* = ui(x) and that 

tr((u(x)u(x)^)“^) = tr(r;(x)“^) = 2£x + 2tr(C'C'*) 
is also a polynomial and therefore in L^(trtc(x)). Furthermore, we calculate 
u(x)“^a 2 (x) + (u(x)“^)'a 2 (x) + u(x)“^ai(x) + ao(x) = 0 
for ao(x) = B. This completes the proof. □ 

Thus as a consequence of the main theorem v = dvi{x) — vo{x) dehnes a Darboux 
transformation from the Bochner pair (e“* 1,6) to the Bochner pair {w,5), where 


4x2/ + CC* 2xC 


, 5 = —d^I + d xl — 


01 -Ai{c*y 


The Main Theorem also tells us a sequence of orthogonal matrix polynomials for w, and 
this allows us to determine a generating function for them. 

Proposition 4.1.2. The monic orthogonal matrix polynomials p{x,n) for the weight 
matrix w define above are given by the generating function formula 

W (cth \ 

^q\C-^ -2{C*C)-^n J nl \ C-H -2{C*C)-^xt - I J 

Proof. The (monic) Hermite polynomials hn{x)I define a sequence of monic orthogonal 
2 

polynomials for I, and satisfy the generating function formula 

j.n 

—,hnix)I = exp(xt — 14)1. 

n 

n=0 

By the Main Theorem, hn{x)I ■ v dehnes a sequence of orthogonal matrix polynomials 

for the weight matrix w. Since the leading coefficient of hn{x)I is I, it’s easy to see that 

the leading coefficient of /i„(x)/ • v is 

c = i ^ (C*)-i \ 

\ C-i -2{C*C)-^n ) ' 
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Thus p{x, n) := c~^hnl ■ v is the monic sequence of orthogonal matrix polynomials for w. 
Multiplying both sides of the generating function formula for the Hermite polynomials 
by u and using this expression for p{x, n) then results in the generating function formula 
in the statement of the theorem. □ 


For N = 2 and C = c ^ C\{0}, this weight matrix appears first in [8| and later 
in [^, where explicit generators and relations of the associated alg ebra D{w) are listed 
without proof. This set of generators and relations is verihed in 28| in a 30-page tour-de¬ 
force. The center of D{w) is also determined explicitly, though misidentihed as being an 
elliptic curve, rather than a singular cubic plane curve. In the following, we demonstrate 
the utility of the main theorem by rederiving the structure of D{w) and the center of 
D{w) succinctly. Better yet, we show that D{w) is naturally identihed with a certain 
subalgebra of 2 x 2 matrices over a polynomial ring, as shown in Equation ([7]). The Main 
Theorem combined with Proposition (11.0.51) gives us a means to calculate the structure 
of the algebra D[w) associated with the weight matrix w. We will use this to calculate 
the structure of the algebra D{w) in the case that N = 2 and C = c ^ C\{0}. 


Proposition 4.1.3. Suppose N = 2 and C = c £ C\{0}. The algebra D[w) is given by 


D{w) 



/ii(e) /12(e) 
/21(e) /22(e) 


u : fij{z) € C[z], 


/i2(|cp+2)=0, /2 i(|c|2)=0, 

/ii(|c|2+2)=/22(|c|2) 


( 7 ) 


for e = — 2dx. 

Proof. Note that D{w) = Mjv(C[e]) for e = — 2dx the Hermite operator. Thus 

D{w) D {iy~^f{e)iy : fie) € MAr(C[e]), ker(i/) • /(e) C ker(i 2 )}. 

If if I and 'i /’2 form a basis for the eigenfunctions of e with eigenvalue |6p, then and 
' 1^2 form a basis for the eigenfunctions of e with eigenvalue |6p -|- 2, and 


i,= ( 

\ A 2 bfj2 J 

is a generator for the kernel of n. Therefore 

D{w) fy Diw, u, w) = {wV(e)i" : /(e) € MAr(C[e]), ■ f{e) G MAr(C)i/)}. 

Now consider /(e) = (fijie)) € M 2 (C[e]). An elementary argument using the fact that 
■01, ) V’l; ^^2 leLust all be linearly independent, shows us that the only way that can be 

an eigenfunction of /(e) is if fuHb]"^ + 2) = 0, / 2 i(|^P) = 0, and /ii(|6p -^2) = / 22 (|&P)- 
Therefore 


Diw, V, w) 



/11(e) /12(e) 

/21(e) /22(e) 


V : fijiz) eC[z], 


fi2{\b\H2)=0, /2i(|b|2)=0, 
/ll(|6|2+2)=/22(|bP) 
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To finish our proof, we must show that the inclusion D{w,i',w) C D{w) is in fact 
an equality. Since the leading coefficient of v is nonsingular, if ry E is such that 

is a differential operator, then the order of is the same as the order of ij. 

( /11(e) /12(e) 

^ /21(e) /22(e) 

s, we see that 


Furthermore the order of 


is equal to 2 maxjj deg(/jj). Looking at 


our conditions for the /, 




D{w,v,w) = Cl ® V I ( ^ ^ 1 /© 1 / ^M 2 (C[e]p(e)p(e — 2))z^, 

V p(e)C p(e)C ) 

for p(e) = e — |6p. This shows in particular that D{w,v,w) has a one-dimensional 
subspace of operators of order 0, no operators of odd order, and a four-dimensional 
subspace of operators of order 2m for every integer m > 0. In other words 


/ /~ \ I 4, if z > 0 is even 

d\m{D{w,v,w)i/D{w,i',w)i-i) = < 

[0, if z > 0 is odd 

Therefore by Proposition 11.0.51 D{w) is generated over D{w, v, w) by operators of order 
0. However, the only operators of order 0 in D{w) are those in C/, and therefore 
D{w) = D{w,v,w). □ 

The proof of the previous proposition shows us that D{w) is given by 


~ 1 ( pi^ ~ 2)C p(e — 2)C \ 1 

D{w) = CI®i' I ^M 2 (C[e]p(e)p(e — 2))i/, 

V p{e)C p{e)C J 

A quick check shows that this algebra is generated by the four elements of order two. 
Moreover, since the center of M 2 (C[e]p(e)p(e — 2)) is C[e]p(e)p(e — 2)1, we have that the 
center of D{w) is 

Z{w) = Cl + v~^C[e]p{e)p{e — 2)Ii>. 

This is a commutative algebra whose spectrum is isomorphic to the nodal cubic curve 
— 2(|5p + l)xy + |6p(|6p + 2)x^ — = 0, via the isomorphism induced by the ring 

isomorphism defined by 

C[x, y] —Cl + i'~^C[e]p{e)p{e — 2)Iv, x 1 -^ v~^p{e)p{e — 2)Iv, y 1 -^ i'~^ep{e)p{e — 2)Iv. 

The ring D[w) has the structure of a module over its center, and as a module it is finitely 
generated and torsion-free. However, it is not a free module over its center since 


0 0 
y(e) 0 


(p(e)p(e - 2 ) 2 /) + 


0 0 
p{e)p{e -2) 0 


(-p(e)y(e - 2 ) 1 ) = 0 /. 


This completes our example. 
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4.2 A Family of Examples of Jacobi Type 

Next, we consider Darboux transformations of Bochner pairs {w, <5) for 
S = —d^{l — x^)I + d{r + 2)xl + B, w = {1 — 

where here 

~ \ 01 T2(/ + 52) ) ' 

for some nonsingular, Hermitian S,T G Mjv/ 2 (C) satisfying = rl, and for 

r > 2. Note that for this to work, we must again assume N is even. According to the 
Main Theorem, we can find examples by finding an invertible, degree 1 polynomial v{x) 
satisfying a certain differential equation, along with a handful of other properties. Define 
ui(x) and uo(x) by 


ui(x) ^ 


s 


0/ 

-S 

01 

T i \ 

-T 

01 


(1-x^) \ vo{x) 


S 0 \ 

0 T J ' 


Then ui(x), uo(x) satisfy the requirements stated in Theorem 1.0.3. for the Bochner pair 

(w,6) 

Proposition 4.2.1. The functions uo(x),ui(x) as defined above satisfy the assumptions 
of the Main Theorem for the Bochner pair {w, S) also defined above. 

Proof. We calculate 


ui(x) 


S-^ 

01 \ ( 

0/ 

-T-i \ 





01 

-T-i i I 

0/ j 


In particular Vi{x) is a polynomial of degree i for i = 0,1 and det(ui) is nonzero on 
(—1,1). We also see that (ui(x)“^)f = (ui(x)“^)* and that 

tr((ui(x)a 2 (x)i;i(x)^)“^) = tr(xi(x)“^a 2 (x)(uF^)*) = tr(S'5* + TT*) ^ 

1 — x^ 

which is in G Lp‘{iT{w{x))dx) since r > 2. Moreover, v^vf^ is Hermitian and therefore 
= VQvf^. Lastly, we calculate 

u(x)“^a 2 (x) + (u(x)“^)'a 2 (x) + u(x)“^ai(x) + ao(x) = 0 

for ao(x) = B. This completes the proof. □ 

Thus as a consequence of the main theorem v = dvi{x) — vo{x) defines a Darboux 
transformation from the Bochner pair {w,6) to the Bochner pair {w,6), where 


51 = (1 - x2)’’/2-i 


S^x^ + T^ -{S^ + T^)x\ 

-{S^ + T^)x T^x^ + S^ ) 
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The Main Theorem also tells us a sequence of orthogonal matrix polynomials for w, and 
this allows us to determine a generating function for them, consistent with the generating 
function found in for N = 2. 

Proposition 4.2.2. A sequence of (non-monic) orthogonal matrix polynomials p{x,n) 
for the weight matrix w define above are given by the generating function formula 

OO 

'^Pn{x)t'^ = 'lpx{x,t)vi{x) + 'lf{x,t)vo{x), 

n =0 

where t) = 2^~^4>{x,t)~^{{l + 4>ix,t))‘^ — and 4>{x,t) = (1 — 2xt + 

Proof. The Gegenbauer polynomials hn{x)I define a sequence of monic orthogonal poly¬ 
nomials for w{x), and satisfy the generating function formula 


'^hn{x)n = 'if{x,t)I. 

n=0 


By the Main Theorem, hn{x)I ■ v defines a sequence of orthogonal matrix polynomials 
for the weight matrix w. Multiplying both sides of the generating function formula for 
the Gegenbauer polynomials by v leads to the desired generating function formula. □ 


For iV = 2, = p, and = r — p, this weight matrix appears in ^ 31|, where 

explicit generators and relations of the associated algebra D{w) are computed, though 
the effort involved is again substantial. In this example, we verify the calculation of 
the algebra D{w) calculated in the paper, using the framework established above. As 
with Tirao’s example (eg. the example of Section 14.11 with N = 2), the machinery gives 
the algebra D{w) with significantly less wrangling. In particular, the Main Theorem 
combined with Proposition (|1.0.5I) gives us a means to calculate the structure of the 
algebra D{w) associated with the weight matrix w, which we do in this special case. 


Proposition 4.2.3. Suppose N = 2 and = p and = r — p for p € (0, r). The 
algebra D{w) is given by 


D{w) = D{w, i',w) = n 


/ii(e) /12(e) 
/21(e) /22(e) 


. fl2(p{n-p+l))=0, f2ii{n-p){p+l))=0, 
fll{p{n—p+i-))=f22{{n-p){p+l)) 


V. 


for e = —5^(1 — x^) + dx{r + 2). 

Proof. Let e = —9^(1 — x^) -|- dx{r + 2). Then D{w) = M 2 (C[e]). Let V’l) V ’2 be eigen¬ 
functions of e with eigenvalues p{r — p + 1) and (r — p){p + 1), respectively. Then 
4>2 = ^n-p '^1 eigenfunciton of e with eigenvalue (r — p){p + 1), and 

(pi = -|- ^-^x'ip 2 is an eigenfunction of e with eigenvalue p{r — p -|- 1). Moreover 

the kernel of u is 


ker(z/) = M 2 {C)'ip{x), ^^{x) 


V’l V’2 
V’l V’2- 
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Our criterion for v to be a differential operator then implies that 

vD{w, u, = {/(e) G M2(C[e]) : V' • /(e) G M2(C)V’}. 

The functions V'l, V’ 2 , ,'/2 are linearly independent. Suppose that 


fie) 


/11(e) /12(e) 
/21(e) /22(e) 


€ M2(C[e]) 


satishes 


"0 • /(e) = Ai/, A = 


a b 
c d 


G M2(C). 


Using the linear independence of the eigenfunctions in the expression 1 / • /(e) = X'ljj, we 
find that 


fi2ip{r-p+l)) = 0, f2iiir - p)ip + 1 )) = 0, fu{p{r-p+l)) = f22iir - p)ip + 1 )). 
From this we see that 

ni~ 1 -^1 ( 

D{w,i2,w) = u { 

[ V /2UO 

Written another way, 

D{w,iy,w) = Cl ® i ^ ^ I 0 z^“^M2(gi(e)g2(e)C[e])z^. 

V 92 (e) 92(e) J 

for g'i(e) = e — p{r — p + 1) and g 2 (e) = e — (r — p)(p 0 1). Since the leading coefficient 
of ly is nonsingular, conjugation by 12 is order-preserving. Furthermore, the order of the 
2x2 matrix (/ij(e) is 2maxjj deg(/ij). Therefore we see that 


/12(e) 

/22(e) 


. /l2(p(?--p+l))=0, /2l((?'-p)(p+l))=0, 
/ll(pT-P+l))=/22((r--p)(p+l)) 


dim(T)(r(;, ly, w)i/D{w, ly, w)i-i) 


4, if i > 0 is even 
0, if z > 0 is odd 


One may check by hand that there are no ui-symmetric operators of order 1, and that the 
only operators of order 0 in D{w) are Cl. Thus by Proposition 11.0.51 we have equality 
D{w) = D{w,iy,w). □ 

It is evident from the presentation of the algebra D{w) given that D{w) is in fact 
isomorphic to the algebra found in our Hermite example (see Proposition 14.1.81) . In 
particular the center is a nodal cubic degenerating to a cuspidal cubic as p ^ n/ 2 . 

The proof of the previous proposition shows us that D{w) is given by 


D{w) = CI + iy \iy + iy ^M2{C[e]p{e)p{e - 2 ))iy, 

V P(e)C p(e)C y 
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A quick check shows that this algebra is generated by the four elements of order two. 
Moreover, since the center of M 2 (C[e]p(e)p(e — 2)) is C[e]p(e)p(e — 2)1, we have that the 
center of D{w) is 

Z{w) = C/ + v~^<C[t]p{e)p{e — 2)Ii>. 

This is a commutative algebra whose spectrum is isomorphic to the nodal cubic curve 
— 2(|5p + l)xy + + 2)x^ — = 0, via the isomorphism induced by the ring 

isomorphism defined by 

C[x, y] —)• CJ + i'~^C[e]p{e)p{e — 2)Iv, x v~^p{e)p{t — 2)Iv, y i-;- v~^tp{e)p{e — 2)11/. 

The ring D{w) has the structure of a module over its center, and as a module it is finitely 
generated and torsion-free. However, it is not a free module over its center since 


0 0 
y(e) 0 


(p(e)p(e - 2)2/) 


0 0 
p(e)y(e - 2) 0 


(-p(e)p(e - 2 ) 1 ) = 01 . 


This completes our example. 


5 General Structure Results 


For the final section of the paper, we will introduce some results regarding the algebraic 
structure of D{w) under fairly general assumptions. To prove the desired results, we 
will use the fact that D{w) is closed under an anti-involution f and that D{w) may be 
embedded into a matrix algebra and is therefore a Pl-ring. This latter fact will allow us 
to apply the results of and 24 ]. 


Before proving our main result, we establish a very helpful lemma. 


Lemma 5.0.4. The algebra D{w) is a semiprime PI-algebra. The nonzero, w-symmetric 
elements of D{w) are not nilpotent. 

Proof. Let {p{x,n)} be a sequence of orthogonal polynomials for w. We claim that if 
6 € D{w) is a nonzero element, then is also nonzero. To see this, suppose that 5^0. 
Then there exists an integer n > 0 such that p{x, n) ■ 5 ^ 0. It follows that 

0 7 ^ (.Pn ' ^jPn ‘ ^)w — (Pn ' j Pn) ■ 


Hence we have that pn ■ 5(5^ is nonzero, and therefore (55^ ^ 0. 

Next suppose that y € D{w) is a nonzero, nilpotent, tc-symmetric element. Then 
there exists a least integer m > 0 satisfying y™ = 0. Clearly m > 1. If m is even, then 
we may write 

0 = y™ = y™/2^m/2 ^ 


By the result of the previous paragraph, this means y”*/^ = 0, and since 0 < m/2 < 
m this contradicts the minimality of m. Therefore m must be odd. However then 
^m+i _ g therefore by the same argument = 0. Since m > 1 we have that 
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0 < {m + l)/2 < m so this again contradicts the minimality of m. Since m is neither 
even nor odd, this is a contradiction. We conclude that nonzero, tc-symmetric elements 
of D{w) cannot be nilpotent. 

Lastly, suppose that X is nontrivial nilpotent two-sided ideal of D{w). Since I is 
nontrivial, we may choose 6 € I with 6^0. Therefore € / is also nonzero. However, 
since I is nilpotent 55^ should also be nilpotent, contradicting the result of the previous 
paragraph. Hence D{w) has no nontrivial, nilpotent, two-sided ideals. This shows that 
D{w) is semiprime. Lastly, the fact that the eigenvalue homomorphism embeds D{w) 
in M]\r{C[n]), combined with the fact that a matrix algebra is a Pl-algebra shows that 
D{w) is a Pl-algebra. □ 

We are now ready to state and prove our general result on the structure of D{w). 

Theorem 5.0.5. Suppose that D{w) contains a differential operator of positive order 
with nonsingular leading coefficient. Then the algebra D{w) is finitely generated as a 
module over its center Z{w) and Z{w) is a reduced algebra of Krull dimension 1. 

Proof. The Krull dimension of D{w) is bounded by the GK-dimension of D{w). Con¬ 
sider the image E{w) of D{w) under the eigenvalue homomorphism S : D(w) = E{w) C 
Misf(C[n]). The Krull dimension of E(w) is bounded by the GK-dimension of E(w) as a 
graded vector space, graded by degree in the variable n. Since E{w) is a graded subalge¬ 
bra of M^r(C[n]), the GK-dimension of E(w) is further bounded by the GK-dimension 
of Mjv(C[n]). However, this latter algebra is Morita equivalent to C[re] and therefore has 
the same GK-dimension as C[n]. Since the GK-dimension and Krull dimension of a com¬ 
mutative ring agree, this means that the GK-dimension of Mi\f{C[n]) is one. Therefore 
the Krull dimension of D{w) is at most one. Since D{w) contains a differential operator 
of order at least one, it must contain a tc-symmetric differential operator 6 of positive 
order d with nonsingular leading coefficient. Then (5"' has order nd for each integer n > 0, 
and it follows that 6 is transcendental over C. Hence the Krull dimension of D{w) is 
at least 1. We conclude that the Krull dimension of D{w) is exactly 1. Combining this 
with the fact that D{w) is a semi-prime Pl-algebra, the main result of tells us that 
D{w) is finitely generated as a module over its center. Since D{w) has Krull dimension 
1, it follows that Z{w) also has Krull dimension 1. 

To show that Z{w) is reduced, suppose that 6 € Z(w). Then € Z(w). To see this, 
suppose that rj G Z{w). Then conjugating everything we find S^rj = rjdfi Since 

r/ G D{w) was arbitrary, this shows that rj G Z{w). It follows that G Z{w), and by 
our previous lemma is not nilpotent. Since it follows that 5™ is 

also not nilpotent. Since S G Z(w) was arbitrary, this shows that Z(w) is reduced. □ 

In the specific case that N = 2, we can say even more about the center Z{w) of 
D{w), namely that it is rational. 

Theorem 5.0.6. Suppose that D{w) contains a differential operator of positive order 
with nonsingular leading coefficient. If D{w) is noncommutative and N = 2, then the 
center Z{w) of D{w) has a spectrum isomorphic to a rational curve. 
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Proof. The previous theorem tells us that the spectrum of Z(w) is a reduced (affine) 
curve. To prove that Z(w) is rational, Liiroth’s Theorem tells us that it suffices to hnd a 
ring monomorphism Z(w) ^ C[n]. We will do so by using the eigenvalue homomorphism 
E : D{w) = E{w) C Mj\f{C[n]). Since D{w) ^ Z{w) we may choose operators 81,82 G 
D{w) such that 8182 7 ^ <^ 2^1 • Moreover we may choose (5i and 82 to be rc-symmetric. 
Then 0 := 1(8182 — 8281 ) is a nonzero rc-symmetric differential operator and hence cannot 
be nilpotent. Setting di{n) := A(5j) and t(n) := A(0), we see that t(n) = i(di(n)d 2 (n) — 
d 2 (n)di(n)) is a trace-free, non-nilpotent matrix in Mf^(C[n]). Since t(n) is trace-free, 
t(n)^ = det(f(n)). Thus since t(n) is not nilpotent, det(t(n))) is not identically zero. 
Thus for all but finitely many values of n di(n)d 2 (n) — d 2 (n)di(n) is a nonsingular 
matrix. Therefore di(n),d 2 {n) generates the full matrix ring M 2 (C) for all but finitely 
many values of n by [l|]. It follows that anything in A(Z(w)) must commute with all of 
M 2 (C), and therefore that A(Z(w)) C C[n]I. Hence Z(w) has rational spectrum. □ 

6 Conclusions 

The goal of this paper is to provide some insight into the rich algebraic structure of 
algebras of matrix differential operators associated to weight matrices. We have demon¬ 
strated a method of constructing numerous examples of solutions to Bochner’s problem 
for matrix differential operators from known solutions. Using our method, we have 
also demonstrated how one may explicitly calculate the associated sequence of monic 
orthogonal polynomials, as well as the corresponding algebra of differential operators. 

We’ve also included some general results regarding the algebraic structure of D(w). 
There are many interesting unanswered questions in this vein that an algebraically in¬ 
clined individual could explore. Examples include 

• Is D(w) equal to the centrilizer in Mjv(U) of one of its elements? 

• Is D(w) isomorphic to the endomorphism ring of a module over its center? 

• Is D(w) maximal among finitely generated algebras over Z(w) with center Z(w)l 

• When is the center of D(w) rational? 

• Is there a stacky picture we can associate to D(w), akin to what we do for heredi¬ 
tary orders? 

We hope to address some of these questions in future work. 
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